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Abstract 

Three fast generators of direct photons in the central rapidity region of high-energy 
heavy-ion colhsions have been presented The generator of prompt photons is based 
on a tabulation of p + p{p) data and binary scaling. Two generators of thermal 
direct photons, for hot hadron gas (HHG) and quark-gluon plasma (QGP) scenarios, 
assume the 1+1 Bjorken hydrodynamics. SPS and RHIC data can be fitted better by 
scenario with QGP. Predictions for the LHC energy have been made. The generators 
have been realized as macros for the ROOT analysis package. 



1 Introduction 

Direct photons are photons not from particle decays. On the quark-gluon level three sub- 
processes dominate: Compton scattering gq 75, annihilation qq 'jg and bremsstrahlung 
emission qq{g) QQig)'!- Photons from initial hard NN collisions are named prompt photons 
and are the main source at large pt- They can be described by perturbative QCD (pQCD). 
In the case of other limit - thermalized system of quarks and gluons, the quark-gluon plasma 
(QGP), these photons are named thermal photons from QGP. 

On the hadron level there are a lot of meson scattering channels: tttt p7, up 
7i'~f,TTK — > K*'-f,Kp K'-f, ... . First two channels give most contribution. If hadron 
system is thermalized these photons are named thermal photons from the hot hadron gas 
(HHG). 

Besides the meson rescatterings, photons from decays of short-living resonances, u 717, 
p — > TTTT'-f, ai Tc-f, A — s> N'j..., make a contribution to direct photons. In the case when the 
life time of a resonance is less than characteristic time of the nucleus-nucleus collision it is 
difficult to reconstruct the resonance because the decay hadron (e.g. tt) can reinteract with 
surrounding medium especially if this medium is dense. 

Below three fast generators of direct photons are proposed: a generator of prompt photons 
and two generators of thermal direct photons. 

2 Generator of prompt photons 

In the paper [1] all existing p + p{p) data on prompt photons in central rapidity region have 
been presented as the function {y/sYEcPa^P/cfp = F(xt),xt = 2pt/y/s. Using a tabulation of 
this dimensionless function one can estimate the prompt photon spectrum in nucleus-nucleus, 
A+B, colhsions at the impact parameter b: Ed^N'^^{h)/S'p = Ed^a^P/d^p ■ AB ■ TAsib), 
where the nuclear overlapping function is defined as TAsip) = Ncouip) / a^^^ where Ncouip) is 
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the average number of binary NN collisions. Nuclear effects (Cronin, shadowing) are ignored 
in this approach. 

Fig. [T] demonstrates reasonable agreement of existing RHIC data [2J with results of the 
generator GePP.C realized as a macro for the ROOT analysis package [3|. Predictions for 
prompt photons at LHC energy is presented in Fig. [2l 

3 Generators of thermal direct photons 

For fast generators the Bjorken hydrodynamics(BHD) [5] has been used. It is assumed that 
during the ion-ion collision the system is mainly expanding in beam direction in a boost- 
invariant way. Natural variables are the proper time r = \/t^ — and rapidity. Thermody- 
namical variables (pressure, temperature, ...) do not depend on the rapidity but are functions 
of r. Viscosity and conductivity effects are neglected. 

Main parameters are the initial thermalization time tq and temperature Tq. Photon yield 
is proportional to ~ Tq, the spectrum slope is controlled by Tq. There is third parameter, 
temperature at freeze-out, Tj. Results depend weekly on it value. Tj=100 MeV has been 
used. 

For this simple space-time evolution one can evaluate expression for photon spectrum with 
the photon emission rate as input which are different for HHG and QGP scenarios. All needed 
formulas can be found in [5]. 

3.1 Generator for the HHG scenario 

The photon emission rates for processes tttt p7, vrp — >■ 717 and p — >■ 71717 the HHG 
scenario have been obtained in [6j using the effective chiral Lagrangian theory with vr, p and 
a\ mesons. Later they were parametrized by formulas in |7j. Rates for the channel a; — > 7r7 
have been taken from [8j. 

Figs.[3]and Hlshow comparison of SPS [9] and RHIC [10] data with results of the generator 
GeTP-HHG.C realized as a macro. SPS data can be fitted by the BHD with the HHG 
scenario, except low points. For fitting the RHIC data unreasonable values of parameters 
should be used, too large ro=7 fm/c, and small To=200 MeV. The HHG scenario is hardly 
realized at the RHIC energy. At the RHIC energy prompt photons dominate science pt=2 
GeV/c while at SPS they are not important up to 4 GeV/c. 

3.2 Generator for the QGP scenario 

In the QGP scenario of high-energy heavy-ion collisions three phases are assumed: a pure 
QGP phase, a mixed phase (QGP and HHG coexist) and a pure HHG phase. Ideal baryon 
free massless parton gas approximation for QGP and ideal massless meson gas approximation 
for HHG are explored. A first order phase transition at the critical temperature Tc is assumed. 
Other parameters of the QGP scenario are the number of degree of freedom in QGP defined 
by the number of colors, A^^c, and fiavors, A^j, and effective number of degree of freedom 
gh in HHG. Time moments between the pure QGP and the mixed phase, r^, between the 
mixed phase and the pure HHG, r^, and time of the freeze-out rj are expressed through the 
parameters. 

The photon emission rates for QGP have been evaluated in perturbative thermal QCD 
applying hard thermal loop (HTL) resummation. Contribution of the next to leading order 
diagrams (bremsstrahlung and annihilation with scattering) is the same order in as from the 
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Figure 1: Direct photon spectrum in central Au+Au collisions at RHIC 
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Figure 2: Prompt photon spectrum in central Pb+Pb collisions at LHC 
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gure 3: Direct photon spectrum for central Pb+Pb collisions at SPS 
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ure 4: Direct photon spectrum for central Au+Au collisions at RHIC 
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Table 1: Characteristics of central collisions at SPS, RHIC and LHC energies in the QGP 
scenario 

leading order diagrams (Compton scattering and qq annihilation). This means that thermal 
photon production in QGP is a non-perturbative mechanism that can not be accessed in 
perturbative HLT resummed thermal field theory. One must consider the obtained QGP 
rates as an educated guess. 

Figs. and M show comparison of the SPS and RHIC data with results of the generator 
GeTP-QGP.C realized as a macro (parameter values Nc=3, Nf=3, Tc=170 MeV and gh=8 
have been used). With the QGP scenario one can try to describe also low pt SPS data. 
The RHIC data can be reproduced with reasonable values of tq and Tq. QGP outshines HHG 
at Pt >2 GeV/c. One can compare the BHD spectra with 2+1 hydrodynamics results |11] . 
While spectrum from QGP is almost the same the HHG spectrum is steeper. One if the 
reasons is the radial flow in the 2+1 case. 

Fig. [7] demonstrates predictions for the LHC energy. Parameter values ro=0.1 and To=650 
MeV have been chosen the same as in 2+1 hydrodynamics [12]. QGP outshines HHG at pt >2 
GeV/c while in the 2+1 hydrodynamics it happens at 3 GeV/c. 

Table [T] summarizes results for SPS, RHIC and LHC energies. Except the initial and time 
evolution parameters the rapidity densities of thermal direct photons, dNj/dy, are pointed 
out. Though QGP dominates at high pt its contribution in dN^/dy is 10%. In the last 
column of the table typical CPU time needed to calculate dN/dpt photon distribution used 
to randomly chose a pt value. 

4 Summary 

Three fast generators of direct photons in the central rapidity region of high-energy heavy-ion 
collisions have been presented The generator of prompt photons is based on a tabulation 
of p + p{p) data and binary scaling. Two generators of thermal direct photons, for HHG 
and QGP scenarios, assume the 1+1 Bjorken hydrodynamics. SPS and RHIC data can be 
fitted better by scenario with QGP. Predictions for the LHC energy have been made. The 
generators have been realized as macros for the ROOT analysis package. First two of them 
have been implemented into the FASTMC code of the UHKM package |13] . 
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Figure 5: Direct photon spectrum for central Pb+Pb collisions at SPS 
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Figure 6: Direct photon spectrum for central Au+Au collisions at RHIC 
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Figure 7: Direct photon spectrum for central Pb+Pb collisions at LHC 
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